Evidence for Scale-Scale Correlations in the Cosmic Microwave Background Radiation 
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We perform a discrete wavelet analysis of the COBE-DMR 4yr sky maps and find a significant 
scale-scale correlation on angular scales from about 11 to 22 degrees, only in the DMR face centered 
on the North Galactic Pole. This non-Gaussian signature does not arise either from the known 
foregrounds or the correlated noise maps, nor is it consistent with upper limits on the residual 
systematic errors in the DMR maps. Either the scale-scale correlations are caused by an unknown 
foreground contaminate or systematic errors on angular scales as large as 22 degrees, or the standard 
inflation plus cold dark matter paradigm is ruled out at the > 99% confidence level. 



Most attempts at quantifying the non-Gaussianity in the cosmic microwave background radiation are motivated by 
the belief that non-Gaussianity can distinguish inflationary models of structure formation from topological models. 
While standard inflation predicts a Gaussian distribution of anisotropies B , spontaneous symmetry breaking produces 
topological defects whose networks create non-Gaussian patterns on the microwave background radiation on small 
scales jg). Minute non-Gaussian features can however be generated by gravitational waves |J or by the Rees-Sciama 
j| and Sunyaev-Zeldovich effects. 

It is generally held that cosmic gravitational clustering can be roughly described by three regimes: linear, quasi- 
linear, and fully developed nonlinear clustering. Whilst quasi-linear and non-linear clustering induce non-Gaussian 
distribution functions, if the initial density perturbations are Gaussian, scale-scale correlations and other non-Gaussian 
features of the density field can not be generated during the linear regime. Hence the linear regime is best suited to 
study the primordial non-Gaussian fluctuations. Since the amplitudes of the cosmic temperature fluctuations revealed 
by COBE are as small as AT/T ~ 10~ 5 , the gravitational clustering should remain in the linear regime on scales 
larger than about 30 hr x Mpc and at redshifts higher than 2. Current limits on non-Gaussianity from galaxy surveys 
probe redshifts smaller than about 1 Q. Interestingly, at redshifts between 2 and 3, and scales on the order of 40 
to 80 h^ 1 Mpc, there are positive detections of scale-scale correlations in the distribution of Lya absorption lines in 
quasar spectra ||. These clouds are likely to be pre-collapsed and continuously distributed intergalactic gas clouds, 
and are therefore fair tracers of the cosmic density field, especially on large scales 0. This may indicate that the 
primordial fluctuations are scale-scale correlated. 

While on small angular scales (£ ~ 150) there may be some indications of non-Gaussianity studies by traditional 
non-Gaussian detectors have concluded that there is no evidence of non-Gaussianity in the cosmic temperature 



fluctuations on large scales ||. (See however |10|.) This does not rule out the existence of scale-scale correlations. 
Because each non-Gaussian feature is non-Gaussian in its own way, there is no single statistical indicator for the 
existence of non-Gaussianity in data. For instance, there are models of scale-scale coupling which lead to a density 
field with a Poisson distribution in its one-point distribution function, but that are highly scale-scale correlated [ pd| . 
In this case, all statistics based on the one-point functions will fail to detect the scale-scale correlations, that is, they 
will miss the non-Gaussianity. As yet, the scale-scale correlations of the cosmic temperature fluctuations have not 
been searched for in any available data set. It is the intent of this Letter to probe for the scale-scale correlations in 
the COBE-DMR 4-year sky maps, and, as an example, show that this measure is effective in testing models of the 
initial density perturbations. In contrast with other techniques, such as the bispectrum, fL2| , higher order cumulants 
Jj"3| , Minkowski functionals jl4| , or double Fourier analysis , scale-scale correlations are localized, and can localize 
the areas on the sky where the signal comes from, and with a resolution that depends on the scale considered. 

The scale-scale correlations are conveniently described by the discrete wavelet transform (DWT) ||[l{|- Considering 
a 2-dimensional temperature (or density) field T(x), where x = (xi,X2), such that < x±,X2 < L, the DWT scale- 
space decomposition of the contrast AT(x)/T is 
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where ipj. ,j 2 ;ii,h( x ) Uiifa = 0, 1, 2.. and l\ = 0, 1...2-? 1 — 1, I2 = 0, 1...2 5 ' 1 — 1) are the complete and orthogonal wavelet 
basis [O. The indexes (ji,^) and (hjh) denote the scale (L/2- 71 , L/2^ 2 ) and position (Lii/2 31 , in phase 
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space and J\ and J2 are the smallest scales possible (i.e., one pixel). The wavelet basis function, ipj 1 j 2 -i ly i 2 (x), is 
localized at the phase space point (ji, ji'i h, h) and the wavelet coefficients £j 1 j 2 ;h.i2 measure the 2-D perturbations 
at the phase space point (ji,j2',h,h)- To be specific, we will use the Daubechies 4 wavelet in this paper, although 
the results are not affected by this choice so long as a compactly supported wavelet basis is used. 
To measure correlations between scales and (j + 1, j + 1), we define 
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where p is an even integer, 1 = (hjh), and the [ ]'s denote the integer part of the quantity. Because Ll/2^ = L2l/2 3+1 , 
the position I at scale j is the same as the positions 21 and 21 + 1 at scale j + l. Therefore, Cj' p measures the correlation 
between scales at the same physical point. For Gaussian fields, Cj' p = 1. Cj' p > 1 corresponds to a positive scale-scale 
correlation, and CJ' P < 1 to the negative case. One can also show that a CJ' P > 1 field cannot be produced by a 
Cj' p < 1 distribution in a Gaussian background. 

It is also possible to define the more "standard" non-Gaussian measures with the wavelet coefficients. Namely we 
define the third and fourth order cumulants as 
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and ej,j;i u i a is the ensemble average (simulated samples) or the average over (hyh) (real data). 

The COBE-DMR data is formatted such that the entire sky is projected onto a cube with each of its 6 faces pixelized 
into 2 10 approximately equal-area pixels. Although one could think of performing a spherical wavelet analysis directly 
on the sky, the current format is ideal for a direct 2-D DWT analysis. The pixels of each face can be labeled by (ji, j'2) 
with < ji, j2 < 5, and I2) with < l\ < 2 31 and < I2 < 2 J2 . The scale j corresponds to angular scale 2.8 x 2 s-3 
degrees. In this way, one can analyze each face individually. This is important, as we can reduce the influence of the 
galactic foreground contamination by selecting the faces in the direction of galactic poles. The galactic plane stretches 
across Faces 1 - 4 (in galactic coordinates) of the projected cube, while Faces and 5 are relatively free of galactic 
interference. We will concentrate on these two faces since the standard galactic cut at \b\ = 20 degrees implies that 
the other faces will be significantly contaminated. 

Before attempting to measure the non-Gaussianity in the DMR maps, we should test for possible contamination due 
to various kinds of noise. A typical example of non-Gaussianity caused by noise is Poissonian noise. Fortunately, this 
type of non-Gaussianity can be properly handled by the higher order DWT cumulant spectra To quantify any 
non-Gaussianity due to DMR noise, we generate 1000 realizations of the temperature maps for a typical CDM model 
with parameters Q Q = 1, h — 0.5, and fit, = 0.05 and generate the appropriate sky maps at the DMR resolution [18). 
To these maps, we linearly add noise to each pixel by drawing from a Gaussian distribution with the pixel dependent 
variance given by the two different foreground removal techniques, DCMB and DSMB (see |l!| for details.) 

Previous non-Gaussian studies using a genus method and other statistics, have found the four year DMR data to 
be consistent with a Gaussian field || . The evaluation of the genus at different smoothing angles is similar to the 
DWT scale decomposition which is also based on smoothing on various angular scales and suggests that the DWT 
cumulant spectra should give similar results. The results for Sj and Kj of the COBE-DMR foreground removed 
maps and the CDM model are shown in Fig. 1. Using the 1000 realizations of the CDM model, we construct the 
probability distribution for both Sj and Kj. Fig. 1 gives the most probable values of Sj and Kj for the CDM model 
with the error bars corresponding to the 95% probability of drawing Sj, Kj from the CDM model. Fig. 1 also shows 
that Sj and Kj for the DCMB and DSMB data are safely within the 95% range. Therefore, one can conclude that 
no significant non-Gaussianity can be identified from the third and fourth order cumulants. This result is consistent 
with the genus results. Note that contrary to previous studies, we can study the six faces of the cube separately. Fig. 
1 shows that both Sj and Kj are isotropic with respect to the face and 5. 

We can now proceed to the scale-scale correlations. We list the most probable values of C*J' 2 for the CDM+DSMB 
maps in Table 1. Similar results are obtained for the CDM+DCMB maps. Because the CDM model is a Gaussian 
model, all C 2 ' 2 are about equal to 1 as expected. At any scale j, C 2 ' 2 is about the same for face and 5. Therefore, 
the noise from the two foreground removed DMR maps does not cause significant spurious scale-scale correlations and 
are thus suitable for a scale-scale correlation analysis. At the very least the sample is good for a comparison between 

observed scale-scale correlations with the CDM model. 
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The results for Cj ' for the COBE-DMR foreground removed maps are plotted in Fig. 2 and tabulated in Table 
1. The behavior of cf ,2) is markedly different from Sj, Kj, or the CDM+DSMB results. First, cf ,2) for face 
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cannot be drawn from the CDM model with a probability greater than 99%. Second, Cj is not isotropic, showing 
a difference between faces and 5. The DCMB maps show the same behavior. 

(2 2) 

C{ ' describes the correlation between perturbations on angular scales of ~ 22 and ~ 11 degrees, which corresponds 
to comoving scales larger than about 100 h^ 1 Mpc. Because the wavelets are orthogonal, Cj; 2 ' 2 ^ cannot be changed 
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by adding any abnormal process on angular scales less than ~ 10 degrees. The "surprisingly" large value for C{ ' 
cannot be explained by any non-Gaussian process on small scales. We have shown that the errors of the foreground 
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removed DMR maps cannot contribute to C{ ' . 

We also checked for possible contributions to the non-Gaussianity from systematics by doing a similar analysis 
on the systematic error maps. It is unlikely that the detected non-Gaussianity comes from the systematics since 
the non-Gaussianity is on the order of ~ 10~ 5 K, while the contribution to the anisotropy from the systematics is 
estimated to be on order ~ 10 -6 K 2Q|. The analysis of the combined systematic error maps confirms that C 2 ' 2 is 

solidly in the Gaussian regime, i.e., C 2 ' 2 = 1.247±0.375. Moreover, these angular scales are larger than the resolution 
of the DMR instrument. Therefore, unless there are very local foreground contaminations which are overlooked by 
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the two foreground removal methods, the high value for and the anisotropy in C\ is cosmological. 

To check if there could be large-scale foreground correlations overlooked by the COBE-DMR subtraction technique, 
we performed the same analysis on the dust maps generated by a careful combination of IRAS and DIRBE data [^p . 
Depending on the method used to obtain an averaged value for the color excess E(B-V) on a DMR pixel, the C{ 
values range from 0.572 ±0.747 to 0.630 ±0.617. Although these maps show small-scale structure, when averaged over 
scales larger than the DMR pixels (2.8 degrees) any non-Gaussian fluctuation disappears. In addition we checked the 
possibility that the non-Gaussianity was due to anisotropies in the synchrotron emission by analysing an all-sky map 
at 408 MHz [Bjl . A visual inspection of this map shows a structure extending from the galactic plane on to the North 
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Galactic Pole. However, using a map projected in the same way as the DMR maps we obtain C{ — 0.069, which 
is much less than the value C{ = 1.008 ± 0.342 obtained by 1000 bootstrap random realizations. As mentioned 
above, C* 2 ' 2 > 1 cannot come from a superposition of a distribution with C 2 ' 2 < 1 in a Gaussian background. Thus 
the scale-scale correlation detected in the COBE-DMR data is not a result of this signal. Additionally, none of the 
individual frequency maps nor a linear combination consisting of the 53 GHz and 90 GHz frequencies maps show 
Ct' > 1.5. Since these maps do not contain the foreground subtractions, this result implies that if the cause of the 
signal in face is foreground, it is incoherent. 

As a final check, we also looked at the correlated noise maps in COBE-DMR ]22|] . The individual correlated noise 
maps of the frequencies were checked for scale-scale correlations and once again, C 2 ' 2 was solidly in the Gaussian 
regime with C 2 ' 2 = 1.007 ± 0.440 for the 31 GHz channel, C 2 ' 2 = 0.896 ± 0.308 for the 53 GHz channel, and 
C* 2 - 2 = 0.935 ± 0.310 at 90 GHz. 

If indeed we have eliminated all non-cosmological sources that could account for this signature and if the signal 
is not just a statistical fluke (since there is still a 1% chance of this occuring), then the only conclusion left is that 
the correlation is cosmological in origin. Whether this signature arises from previously proposed sources of non- 
Gaussianity, such as cosmic strings, large spots, matter-antimatter domain interfaces, etc., remains to be determined. 

Recall that the COBE-DMR data tolerate almost all popular models of primordial density perturbations in terms of 
the second order statistics. Generally, the data are only able to discriminate among the power spectra of these models 
with less than 2-a confidence levels fl24[ ]. The scale-scale correlation detected in the 4-year COBE-DMR data either 
gives a rather high confidence of ruling out the CDM model or evidence for the existence of unknown local foreground 
contamination on angular scales as large as ~ 10 - 20 degrees. Obviously, if either of these implications are correct, 
two important conclusions can be inferred: 1.) the inflation plus cold dark matter model with standard cosmological 
parameters appears to be ruled out at the > 99% confidence; 2.) the COBE-DMR temperature maps are contaminated 
on large angular scales at levels larger than previously thought. Whether COBE determined cosmological parameters, 
such as the quadrupole of temperature fluctuations, may also be contaminated remains to be seen. 

We are very grateful to Al Kogut for providing the systematic error and 408 MHz maps used in the analysis. 
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TABLE I. Measured Cf 2 coefficients. 



FACE 


FACE 5 


3 


DMR 


CDM+DSMB 


95% Confidence Levels 


DMR 


CDM+DSMB 


95% Confidence Levels 


1 


2.091 


1.004 


(0.376 - 1.587) 


0.730 


1.008 


(0.492 - 1.761) 


2 


0.984 


1.035 


(0.601 - 1.513) 


1.300 


1.022 


(0.688 - 1.572) 


3 


1.041 


1.032 


(0.791 - 1.294) 


1.172 


1.026 


(0.832 - 1.358) 



4 



22° 11° 5.5° 



to 



- 


f 


I 


* 




<> 


,* 




£ 


I 

) 


-A 


- 




A 




! ▲ 




A 




















Face 


: 


o CDM 


A 


DCMB 


£ DSMB 










Face 


5 : 


CDM 


▲ 


DCMB 


★ DSMB 
















[ 




c 


) 


▲ 

★ 






c 

A 


> 


★ 




•k 

A 


▲ 

* 




* 












I 














I 







2 3 



FIG. 1. Sj (top) and Kj (bottom) for faces and 5 of the DMR data and the CDM simulations. 
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FIG. 2. C (2,2) for faces and 5 of the DMR data and the CDM simulations. 



